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Abstract. We use novel integral representations developed by the sec- 
ond author to prove certain rigorous results concerning elliptic boundary 
value problems in convex polygons. Central to this approach is the so- 
called global relation, which is a non-local equation in the Fourier space 
that relates the known boundary data to the unknown boundary values. 
Assuming that the global relation is satisfied in the weakest possible 
sense, i.e. in a distributional sense, we prove there exist solutions to 
Dirichlet, Neumann and Robin boundary value problems with distri- 
butional boundary data. We also show that the analysis of the global 
relation characterises in a straightforward manner the possible existence 
of both integrable and non-integrable corner-singularities. 



1. Introduction 

In this paper we present a new means to study rigorous aspects of bound- 
ary value problems for elliptic PDEs in convex polygons. We concentrate 
on the basic elliptic equation 

(1) - Ag + fi\ = 0, 

where > (modified Helmholtz). The case /3 = (Laplace) was studied 
in [9] for continuous boundary data and can be extended to distributional 
boundary data using the methods presented here. All our results hold if 
< (Helmholtz or Laplace), but for economy of presentation we state 
our results for > only. 

We study (a) existence of solutions with distributional boundary values, 
(b) a priori estimates for boundary data regularity and (c) corner singu- 
larities. The approach we use is based on the formal results obtained in 
[6] where novel integral representations for solutions to the basic elliptic 
equations in a convex polygon were derived, under the assumption that the 
solution exists. The rigorous methodology is based on the following pro- 
cedure: we define a function using the novel integral representations of \Q\ 
and show that this function satisfies the relevant PDE and converges to the 
desired data on the boundary provided that the global relation is satisfied. 
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Regarding (a), we first generalise the integral representations in jfi] so 
they make sense for distributional boundary values. We then show that 
the function defined by the associated integral representation satisfies the 
boundary value problem with distributional boundary values provided that 
the unknown boundary values satisfy the global relation in the distributional 
sense. 

Regarding (b), we provide classical a priori estimates for the unknown 
boundary data by analysing the global relation. 

Regarding (c), we use the global relation to prove the existence of corner 
singularities when certain mixed boundary data are prescribed at adjacent 
edges of the polygon. The arguments are based on simple asymptotic eval- 
uation of terms in the global relation. 

In this paper we do not address the fundamental question of existence 
of solution to the global relation. It is known that if Dirichlet data is in 
if*^^/^(5r2) with s G [— ^, then the global relation uniquely defines the 
unknown boundary values, which belongs to [2]. A new, more 

direct method for solving the Dirichlet-Neumann map via the global rela- 
tion was presented in [3] for Dirichlet Data in H^{dQ). This method can 
be extended to boundary data with lower regularity. We expect that if one 
works on a suitable quotient space, then regularity results for the Dirichlet- 
Neumann map will hold for general distributional boundary data. The rele- 
vant space quotients out finite linear combinations of Dirac measures at the 
vertices (the relevance of this space is discussed in Theorem [T]) . 

For a more classical approach to boundary value problems with distribu- 
tional boundary data, we refer the reader to [H [3] and references therein. 

Several numerical techniques for the implementation of the unified method 
of [U [7] applied to linear elliptic PDEs are discussed in [HI |9l [151 ES] ■ 

2. The Spaces of distributions 

Let Q, C ~ C denote the interior of a polygon with vertices {zi}f^i 
and sides Fj = {zi,Zi+i). Let Oj = arg(zj4.i — Zj) denote the angle the side 
Fj makes with the positive real axis and let |Fj| be the length of the side 
Fj. It was shown in [6] that if one assumes that there exists a solution 
to ([T]) in J7 with sufficient smoothness, then the solution has the following 
representation: 

(2) ^(-'-") = 4^E/^-^'^'^'^'/V.(A)^, 

where the functions {pi{X)}^^i, called spectral functions, are defined by 



^'l \ A ' ^ f^^^ ^'i \a-' 
- + iAg d. + -,--)dz 



and {fj}"^]^ denote the rays in the complex plane with arg(A|^-) = —a 



(3) pi{X) = / e-'^-'+'^ ^"'/A 

and {fj}"^]^ denote the rays in 
orientated towards infinity. Furthermore, the spectral functions satisfy the 
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global relation 

n 

(4) Y.p,(X) = o. 

1=1 

In what follows we will first generalise these results so that they make sense 
for distributions (generalised functions) . In this respect we will use standard 
results from the theory of distributions, a comprehensive account of which 
can be found in [lOj. We will be interested in distributions supported on the 
closure of the edges Fj (the edges Fj do not contain the vertices). We use 
£"[0, 1] to denote the space of smooth functions from [0, 1] to C. This space 
is endowed with a topology generated from the semi-norms 



Sn{^) = max v?(")(r) 

rG[0,l] 



n = 0,1,2,, 



a sequence of functions (/^^ — >• in £[0, 1] if and only if Sni^Pk) ~^ for all 
n G N U {0}. This is a Frechet space and we denote its topological dual by 
£^'[0, 1], which is the space of continuous linear forms from S[0, 1] to C. The 
pairing between £^'[0, 1] and £^[0, 1] is denoted by (•, •). It is well known that 
a linear form u : £^[0, 1] — )■ C is an element of £'[0, 1] if and only if there are 
constants N and C such that 

(5) \{u,ip)\<CY,U^)^ 

n<N 

for ah (pG£[0,l]. 

We will also make use of the less familiar space of Schwartz distributions 
on an open interval. To define this space of distributions, we first setup the 
appropriate space of test functions. Let 5(0, 1) denote the space of smooth 
functions from (0, 1) to C such that 

amn{^)= sup |r-™(l-r)-V"HT)| <oo, n,m = 0,1, 2,..., 

re(0,l) 

so that these functions vanish to all orders at the boundary of [0, 1]. The 
topology generated by these semi- norms makes 5(0, 1) a Frechet space. The 
topological dual, consisting of continuous linear forms from 5(0, 1) to C, is 
denoted by 5'(0, 1). A linear form v : 5(0, 1) — )• C is an element of 5'(0, 1) 
if and only if there exist constants D and M such that 

m,n<M 

for all if G 5(0, 1). For if E 5(0, 1) we define the Fourier transform by 

^(A) = / e^'^Xr) dr. 
Jo 

If (p{x) = (p{—x), the Fourier inversion theorem can be written as 27r(^ = 
((i^)")". Using the natural extension of an element of £^'[0, 1] to £^'(R), the 
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space of distributions on R with compact support, we can define the Fourier 
transform on £'[0, 1] via 



See |10l Ch. 7] for the standard treatment of the Fourier transform on iS'(R). 

We now generaUse the integral representation ^ so it makes sense for 
distributional boundary values. To facilitate this generalisation we need a 
local description of the edges Fj so that we can define distributions on them. 
We use the local parametrisations ipi : [0, 1] — >• Fj, with 

Ipi'.T^ TZi+l + (1 - T)Zi. 

For a function / : Fj — )• C we write its pullback by ■i/'i via 

V'*(/)(r) = /(V(r)). 
1-forms are pulled back by ipi via 

ipUdz) = e'"'|Fi| dr, ij*{dz) = e-^°»|Fi| dr. 
We now define the spectral functions Pi{X) by 

PiW _ /a „ ,,.*/„-iA2'+i/32j'/A^ 



(6) 



ilFi 



+ ( Ae'"- + 



where dnqi,qi G £'[0, 1] for i = 1, . . . , n. The first term is motivated by the 
identification 

-1 



— dz' - — dz 
dz' dz' 







dq 



ir^l (^(9n%,'0*(e 

In this calculation we have used the identities 



dn 

AXz'+ip'^z' /X 



(r) dr 



dq 
dz' 



dq . dq 

at ^'cT 



n 



dq 



dq . dq 
9t "'cT 



n 



where d/dt and d/dn are derivatives in the tangential and (outward) normal 
directions along Fj. The second term in ^ is similarly motivated, e.g. 

T 



A 

i/32|Fi|e- 



AXz'+ip'^z' /X 



(r) dr 



AXz'+ili^z' /X 



So the formal identifications are 



dnqi ^ ipi 



dq 
dn 
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Figure 1. The boundaries dfl and dfl^. 

3. Distributional Boundary Data 

In this section we will outline how to setup a given boundary value prob- 
lem so that the boundary data can be distributional. Our aim is to show 
that the function defined by ([2]), where the {pi}f^^ are given in ([6]), satisfies 
the following boundary value problem: 

(7a) -A^ + /32g = inn 

(7b) lim q = Qi in S'{Ti) for z = 1, . . . , n, 

where qi £ <f'[0, 1]. We need to specify a notion of convergence, via which 
the solution ([2| converges to a distribution on Fj. 

For a given polygon Q we inscribe within it a one parameter family of 
polygons 

n^ = {z en-. dist{z,dn) > e}. 

For e > sufficiently small we have dQ ~ and we can define the edges 
of dQe by for i = 1, . . . , n as shown in Figure [3j Now suppose q G C°°(f]) 
is a given function. Using a local parametrisation ipf : [0, 1] — F^, we can 
define a distribution g| G 5'(0, 1) by 

(8) {q'i,^)= f m\q){TMT)dT, 

Jo 

for all ip £ 5(0, 1). As e | we have F^ — )• Fj and ipf — )• ipi. This limit 
and its equivalence to others modes of convergence for boundary values 
are discussed in [4J (cf. [H]). Our definition of g E C'^{Q) tending to a 
distribution on Fj is 

I lim g = qi in 5'(Fi)| ^ | limg| = qi in 5'(0, 1)|. 

These limits are understood in the topology of 5'(0,1). Recall that our 
boundary datum is an element of £'[0, 1]; on the other hand we are taking 
limits in the topology of 5'(0, 1) and not in the topology of £'[0, 1] (the test 
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function appearing in ([s]) is an element of 5(0, 1) and not of iS[0, 1]). Thus, 
we must make sure that we do not lose any information in this limit process. 

In this respect we note 5(0, 1) C <S[0, 1] and that the inclusion is con- 
tinuou^ Hence ^'[0, 1] 5'(0, 1). Therefore there is a natural projection 
TT : <?'[0, 1] — >■ 5'(0, 1), which is obtained by restricting u G <f^'[0,l] to test 
functions in 5(0, 1). Given a linear form v G 5'(0, 1), the Hahn-Banach the- 
orem implies that we can extend this linear form to a linear form uq G S'[0, 1] 
with tt{uq) = V. The general solution to the equation tt{u) = ?; is |13l p. 
178] 

N N 

k=0 k=0 

This observation means that the difference between u G £'[0, 1] and ir^u) can 
only be a collection of Dirac measures and derivatives thereof supported at 
the boundary of [0, 1]. Intuitively speaking this is because the test functions 
in 5(0, 1) vanish to all orders at the boundary of [0, 1], so they are unable 
to "detect" a distribution in S'[0, 1] with support contained in the boundary 
of [0, 1]. These observations imply that knowledge of 

lim q = Qi in S'iTi) 

yields 

TV 

(9) ]im q = q,+Y, (a^S^ + b,6^ \ in £'{T,), 

" 3=0 

for some constants {aj,6j}^Q and some integer A^. The possibility of the 
existence of the above additional contributions to the right hand side of ([9| 
does not affect the answer to the questions of existence of a solution to ([T]), 
since we can always consider the new function (distribution) 

n N 
i=l j=0 

for appropriate {cij} which satisfies the relevant PDE in and converges to 
the desired distribution in £"'[0, 1]. 

The following result shows that 5'(0, 1) is, in some sense, the right space 
of distributions to take limits in. 

Theorem 1. Let A[0, 1] denote the space of all finite linear combinations 
of Dirac measures and derivatives thereof supported at and 1 . Given qi S 
£^'[0,1] (i = l,...,n), there is at most one of each d^qi G f [0, 1]/A[0, 1] 
(i = 1, . . . ,n) such that the global relation is satisfied. 

This result implies that a solution of the global relation is unique modulo 
Dirac measures at the vertices of the polygon. By taking limits in 5'(0, 1) 

^Meaning that if a sequence {(pk}k>i tends to zero in 5(0, 1) (i.e. amni^Pk) — ^ for all 
m, n), it also tends to zero in £[0, 1] (i.e. 5n{ipk) for all n). 
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we effectively quotient out this degeneracy. In this sense, the topology on 
5'(0, 1) is the right one to take limits in for distributional boundary data. A 
more detailed treatment of distributional solutions to the global relation will 
be presented elsewhere, but the interested reader may look at [3l Lemma 4] 
for an idea of the proof. 

To prove our main result we will need to interchange the order of some 
integrals. To this end, the following simple lemma is essential. 

Lemma 1. For z £ il. and for each i, the function 



A 



decays exponentially as |A| — t- and 



oo along the ray ii 



Proof. From the semi-norm estimates ^ we know that there exist constants 
C, N such that 



P.(A)|<CE(|A| + ^ S 

n=0 11/ 



-iXz'+i/S^z'/X 



The proof now follows from calculations analogous to those in O Lemma 
1]. ■ 

From this observation it is clear that the double integrals 



Jo Je. 



dX 

T 



V) dr 



converge absolutely for each 1 < i, j < n and e > 0. Hence applying Fubini's 
theorem we can interchange the order of integration 



e\* / giAz— i/3^z/A 



^p{t) dr 



dA 



Lemma 2. Let £i denote the ray for which arg(AU ) = tt — Oi, so that ii is 
the continuation of the ray ii . Then, for j ^ i: 



[ pjwir 



^*|giA.-i/32.-/A 

"1 



. pjW 



(P{t) dr 



dA 
T 







(P{t) dr 



dA 



Proof. In order to be able to deform contours in the complex plane, we need 
to establish some estimates for the integrand when |A| is either small or 
large (Cauchy's theorem will then provide the result since the integrand is 
analytic in any punctured disc centred at the origin). In this respect, for 
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each m > and some constant Cm we have the following estimates: 
"1 







< Cm |A| + 



(/3(r) dr 



< Cm |A| + 



|A| 
|A| 



max 

ci{r,) 



^*(e'^^-'^'""/^)y.("^)(T)dr 



This estimate follows from repeated integration by parts, using the fact that 
(p and its derivatives vanish to all orders at the boundary dTi. The semi- 
norm estimates used in Lemma [T] imply that for m and Cm as above, we 
have 



< C„. (|A| + g 

Observe that 
(10) 



(P{t) dr 



max max 

2eci(ri)2'eci(rj) 



,i\{z-z')-i/3^iz-z')/X 



exp 



|A| + ^)sin(arg(A(z-z'))) 



If |ai — aj\ = vr, then £j and ii coincide, thus there is nothing to prove. Let 
us suppose < ai — aj < vr. This condition together with convexity implies 

Oi — TT < aj < arg(z — z') < aj + vr. 

Hence, if A satisfies vr — < argA < —cej, we find 

< arg(A(z - z')) < vr. 



It follows now from (10) that the ray £j can be deformed onto ii. The case 
< Q,- — a,- < vr is similar. ■ 



Theorem 2. Define the function Q = Q{z,z) by 



Q{z,z) 



n 

— Y 

4vri ^ 

4 = 1 



JiXz-ip'^z/X 



dX 



with {/Oj(A)}j'L]^ defined in where qi, dnQi G £'[0, 1] for i = 1, . . . ,n. If 
the global relation Q is satisfied, then Q solves the following boundary value 
problem: 



(11a) 
(lib) 



-AQ + l3^Q = inn 

lim Q = Qi in S'{Ti) for i = 1, . . . , n. 

z^Ti 



Proof. It is straightforward to show that Q satisfies the relevant equation 
at any interior point of O, since exponential convergence justifies differenti- 
ating under the integral sign. Thus, we focus on the behaviour of Q at the 
boundary. Without loss of generality we position the vertex Zi at the origin 
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and align the side Fj with the positive real z-axis so that Qj = 0. Our aim 
is to compute the limit 



lim {Q<r, ^) = lim / (^|)*(Q)(r)(^(r) dr. 

We have already established that the double integral converges absolutely, 
so interchanging the order of integration is justified. Thus, our task is to 
compute the limit 

-1 



.7 — -L 



dA 



It is simple to bound each of the integrands by an absolutely integrable 
function (independent of e), employing the same integration by parts type 
estimates used in Lemma [2} So the dominated convergence theorem applies 
and we can pass the limit through the integral signs. Hence, 



lim (Qt, 



J — I J 



dA 
T 



1 
47ri 



-1 



A(A) 



(P{t) dr 



dA 

T 

(P{t) dr 



dA 
T' 



Using Lemma [2] and the global relation Q we find 



dA 

T 



1 

47ri 



ifir) dr 



dA 

T 



1 
Am 



dA 



Hence 
(12) 



1 

47ri 



a(a) 



* / iAz-i/J^z/A 



(/9(r) dr 



dA 



Since = 0, we can parametrise the A-integrals in (12) so that ii is the 
positive real axis and ii is the negative real axis: 



4vri Wuv iQl^) 

£-5-0 



>(|r,|(f -a))p.(a) 



dA 
T 
dA 
T 



V(|r.|(f -a))a(A) 



<^ ir^l A 



P.(-A) 



dA 
T 
dA 
T' 
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In the first integral we make the change of variables k = |rj|(A — i.e. 

where we have chosen the positive square root so that A > 0. Similarly, in 
the second integral we make the change of variables k = — |rj|(A — /3^/A), 
i.e. 

2A = {-k + y^^TW^' 



ird ^ ' " " y ' A v^fc2 + 4|ri|2/32 

Hence 

' ^A^^\{^r.qir{k) 



47ri lim(Q,S(/^) = / ^{-k) 



(13) - / ^{-k) 



' i|ri|(5„(?,)lfe) 




yA;2 + 4|r,|2/32 

We note that the above integrals are well defined for Qi^dnQi G iS'[0, 1], 
since by the semi- norm estimates ([s]), the corresponding Fourier transforms 
are entire functions of polynomial growth on the real line, . We find the 
contributions from the derivatives cancel, and we are left with the expression 

(14) lini {Qt, ^) = ^ r V>i-mik) dk. 

e^O 27r 

Now if E 5(R) and since qi S £'[0, 1] ^ 5'(R), we can apply the Fourier 
inversion theorem on 5'(R) to find 

lim {Qi, = 1^- {Qi, (0)") = {Qt, ■ 

Hence lim2_5.r, Q = qi m S'{Ti). ■ 

Remark 1. By performing a similar calculation it is possible to determine 
the limit of dQ/dn, where n is the normal to F^ ~ Fj. In this case the 
corresponding terms for in (113^ cancel and we find 



lini {{dQ/dn)l ip) = {dnqi, ■ 

Thus, under the assumption that the global relation is satisfied, it is possible 
by using Q, to prove existence for the Dirichlet, the Neumann or the Robin 
boundary value problems. 

Remark 2. Central to this argument was the statement that our domain 
can be placed in the z-plane so that the vertex at z = Zi lies at the origin. 
We certainly do not lose any generality by enforcing this situation. Indeed, 
suppose a priori that one is interested in the behaviour of the boundary 
values at a particular edge of the domain boundary. Then, after this edge 
has been singled out, the domain can be embedded into C ~ R2 in any way 
we find convenient. Alternatively, one can check directly that the integral 
representation ([2]) is invariant under rotations and translations in the z- 
plane. 
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This highlights an important advantage of the unified method introduced 
in [HI [7] : in ah boundary value problems there is an underlying gauge free- 
dom, namely, the domain on which the boundary value problem is defined 
may be embedded into R" a variety of ways. Clearly, the particular choice 
of gauge does not affect the analysis of the given boundary value problem, 
since the solution at a point x S $7 is determined by the boundary data 
prescribed on d^l, and the position of the point x relative to the boundary. 

In the classical approaches (Green's functions, boundary integral equa- 
tions, etc.), the analysis is carried out in the physical space. In this setting, 
gauge freedom gives no particular advantage (the action of a gauge transfor- 
mation simply changes all variables appearing in the equations). However, 
in the unified approach the analysis is done entirely in the spectral space. 
The spectral parameter is decoupled from the physical parameter, and now 
a choice of gauge does affect the analysis of the spectral problem. In partic- 
ular, as we have seen in this section, a choice of gauge allows one to single 
out a particular neighbourhood of the boundary d^l. Since the underlying 
problem is gauge invariant, these local arguments can be extended to all of 

dn. 

From the above discussion, it follows that the "without loss of general- 
ity" statements made in this section can be interpreted as an exploitation 
of gauge freedom: we have made a particular choice of gauge so that the 
analysis of the boundary values on part of our domain is simple, then we 
have extended these results to the entire domain using the underlying gauge 
freedom. 

4. Regularity of Boundary Data 

In this section we use the global relation Q to derive regularity results 
for the unknown Neumann data. On the assumption of a distributional 
solution to the global relation, we use elementary results from the theory 
of distributions to retrieve the classical regularity results for the Dirichlet- 
Neumman map for elliptic problems. The basic result is summarised below 

Theorem 3. Suppose the distributions qi,dnqi G '^'[0,1] (i = 1, . . . ,n) are 
such that the global relation Q is satisfied: 



(15) >.a(A) = 

where 



n 
i=l 



^^^^^ -/5nft,,V*(e-'^^'+'^'^"'/^ 



Suppose that the natural extension ofui belongs to H^(Fi) for some 1 < i < n 
and sen. Then dnQi G -^^loc^O, !)• 



12 



A.C.L. ASHTON AND A.S. FOKAS 




The main idea in the proof is to use the gauge freedom to ahgn the 
particular side we are interested in with the real z-axis, so that all other 
sides of the polygon are contained in the lower half plane, as in Figure [2] 
The term in the global relation corresponding to the side aligned with the 
real z-axis, is closely related to the Fourier transform of the boundary values 
on this side. By proving that the other terms in the global relation decay 
in a suitable sense, the regularity of the unknown boundary values on the 
side aligned with the real z-axis can be deduced from the regularity of the 
known boundary data using standard properties of the Fourier transform and 
Sobolev spaces. The proof of the theorem requires two technical lemmas, 
which show that if we are interested in the regularity of the boundary values 
on the side Fj, the only important term in the global relation is pi. 

Lemma 3. Let he scaled and positioned so that it resides in the lower 
half z-plane and the side Fj coincides with the interval (0,1) on the real 
axis. Then 

(16) p^{X) = {^nq^nX " ^) " f A + ^X - P^/X). 



X , 

Then the following estimate is valid for A > 0.- 

(17) |/j._i(A) + p^{X) + p.+i(A)| <c[x + ^\ e-<^+P""'\ 
for some fixed e > and constants C,M > 0. 



Proof. The equality in (16) follows immediately from the fact that Fj co- 
incides with (0, 1) and = tt. If j ^ ± 1}, then the basic semi-norm 
estimate ([s]) gives 



|p (A)| < C ( A + ) max e^^+^'Z^)'-^ < c ( A + ^ ) g-^^+^'Z^) 
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since \m z' < —e on Tj for some e > 0, using the fact that all the sides lie 
lower half plane. The estimate in (17) follows from the global relation (15) 

Pi_i(A) + pi{X) + Pi+i{X) = - Pj{X) 

and the previous estimate for the terms on the right hand side. I 

Lemma 4. For if G C^(0, 1) we define the operation * : 5(R) xC°°(R^) — ;> 
C7°°(R) 

poo . 2 \ 

(p * p(A) = y (a - /i + ^ j p{ii) dM(/x), 
where dM{p) = (1 + l3'^/p'^)dp. Thenipi^pi e 5(R). 

Proof. First note that since if G C^(0, 1), we have G ^(R) and then it is 
straightforward to show that 

This also means that this map is well defined, because 



/3 



The estimate in (17) and the closure of 5(R) under convolution imply 

(p-k Pi = -ip-k pi_i - (f-k pi^i modcS(R). 

We need to show that the terms on the right hand side of this expression 
are Schwartz functions. The analysis is similar for both terms, so we focus 
attention on ip-k pi+i. Our object of interest is given by 



(18) = lim / 0[X-p+ ^) p,+i(^)e-('^+/^'/^) dM{p), 



where we have inserted a regularisation for later convenience. This is jus- 
tified since pi-^-i{X) is certainly polynomially bounded (by the estimate in 
lemma [3]) and (p € 5(R), so the dominated convergence theorem applies. 
Now one can show that any distribution in £'[0, 1] can be written as a fi- 
nite sum of derivatives of bounded functions (see appendix). Employing 
this result we may assume that there exist bounded functions {gj}jLi with 
supports contained in [0, 1] such that 

j=oJo V ''"/ 

and similarly for the term corresponding to qi+i- Hence a generic term 
appearing in />i+i(A) can be represented in the form 

/V(r)e-'^^^'+2+^^'^^»+2/^dr, 
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where / is some bounded function and /3 G Z. Since pi+i{X) is just a finite 
sum of such terms, it is sufficient for our analysis to use this expression 
instead of pj-i-i(A). Applying Fubini's theorem, equation (18) reads 



iflX- H 



§1 



13 / yj^^-jg-iM-r2i+2+i/3^TZi+2/At 



dM{fi) 



(19) 
where 



lim 



''^V(y)F,(y)dy, 



fir) 



dr. 



= fir) 



dr, 



with Zeir, y) defined accordingly. We now proceed to show that the function 
F{y) = lime_>.o -Fe(y) is smooth on the interval (0,1). Making a change of 
variables in the //-integral 

p! = flZeiT,y), 

and noting that ReZeiT,y) > e, we can use elementary techniques from 
complex analysis to show that the corresponding contour of integration in 
the //'-plane can be rotated back to the real axis. This gives 
™i /(^) r .oo 



Feiy) 



Zeir,y) 



{p')^e->''e-^\'^^^^y^\/^''dM{p:/z,) 



dx. 



For any y G (0, 1) there exists a (5y > such that \zq{t, y)\ > Sy for r G [0, 1] 
so we can take e — )• for y in this open interval, although the convergence is 
not uniform in y. Also, for any y £ (0, 1) the //' integral decays exponentially 
so we can differentiate under the integral sign as many times as needed. It 
follows that F € C°°(0, 1) and hence cpF £ C~(0, 1). We deduce from ^ 
that ip-kpi^i £ 5(R), since it is the Fourier transform of a smooth function 
with compact support. The argument showing ip-k pi_i G 5(R) follows in a 
similar fashion. ■ 

Proof of Theorem Using the result of the previous lemma we have 

^★/JiG5(R) V(^gC~(0,1). 

We introduce the new variable k = p — 0^ / p for the integral defining <f>-k pi. 
It is clear that k[p) : — t- R is a diffeomorphism and dk = dM{p). Using 
this observation in (|16|) we find 



p^{p{k)) = id^qiYik) - VP+4/32g,(fe). 
Making the substitution in the integral for pi we find 

/oo /"OO 
ip{X-k){dr,qi)Xk)dk- ^p{X-k)y^WTWqiik)dk = mod5(R). 
-oo J — oo 
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The first term on the left hand side it just <p * (dnqi)", so 

{(fdnqiTiX) = ^{^* (dnqdlW 
(20) =— / ^(A- fc)\/^2 + 4^2^.(^)(i^ mod5(R). 

Using the standard notation (A) = (1 + A^)^/^, it fohows that 

{k) ^ {k" + 

We want to estimate the norm of the right hand side of (20) after mul- 
tiphcation by (A)*~^, i.e. we want to estimate the norm of 

—\^{k)\{\-kY\qi{\-k)\dk. 



,{X-k) 
Using Peetre's inequality 



-oo (A- A:) 



and applying Young's inequality we find 

\(p{k)\ (A - kY \qi{\ - k)\ dk < Hy'll/z^-i Ikilk^- 

Using this estimate in (20) we find that ^pdnqi G H^~^(R). Since f G 
C^(0, 1) is arbitrary, we conclude that dnQi G H^^^{0, 1). ■ 

Remark 3. The proof is almost the same if one assumes that if G 5(0, 1). 
The only modification needed is to take into consideration that the function 
F{y) = lim^^Q F^{y) has a singularities of finite order at y = 0. It is then 
easy to show that ipF G 5(0, 1), so its Fourier transform belongs to 5(R) 
and the proof follows as in Theorem |3j 

We note that this proof can easily be modified to deal with more general 
Sobolev spaces VF'^'^(rj) where p G (l,oo). Here we only state the relevant 
definitions and results, and refer the reader to for a comprehensive 
account. 

The class of Fourier multipliers 5'"(R) are defined by the smooth functions 
p{x, A) such that 

|9^a^p(x,A)|<C„;3(Ar-H. 

For p G S''^(R) we define the pseudo-differential operator p{x,D) G \1'"(R) 
via 

poo 

p{x,D)f{x)= / p(x,A)/(A)e'^^dA, 



for Schwartz functions /. The action of p{x, D) can then be defined on 5'(R) 
by duality. A fundamental property of the pseudo-differential operators 
^'"(R) is that for q G (1, oo) the following map is continuous, 

(p, f) G 5"(R) X W''''{-R) ^ p{x, D)f G <r'''(R), 
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meaning that for p G 5"'(R) and / E VF*''^(R) with q G (l,oo), we have 
\\p{x,D)f\\^s-^„<C\\f\\w^.,, 

loc 

for some constant C. We say that D) G ^^"'(R) is elliptic of order n if 

\p{x,\)\>c{\r\ 

for |A| sufficiently large. With these definitions we have the following result. 



Theorem 4. Let {dnqi\^^i and {qi\^^i he distributional solutions of (15). 
Then, there exists an elliptic pseudo-differential operator p{D) G ^""^(R), 
such that for each ip G C^(0, 1) 

(f {dnqi - p{D)qi) =0 mod5(R). 



Proof. Recall the equality ( 20 ) 

■o2 roo 

(21) {v^^q^r{X) = — / k)^/W^Wqt{k) dk modcS(R), 

^ J-oo 

is valid for all if G C^(0, 1). We introduce the function 

p{k) = 2il3^^/WTW- 

It is clear that p{D) G ^^(R) is elliptic. Working modulo the space of 
Schwartz functions, the equivalence in (21) reads 

(<^%)"(A) = ^ r ifiX - k)p{k)q,{k) dk 

= l-0*{p{D)qnX) 
= {MD)q)XX). 

So we certainly have 

[ifdnqi - (pp{D)q,i)XX) = modcS(R). 

The result follows from the fact the Fourier transform is an isomorphism on 
5(R). ■ 

5. Corner Singularities 

We consider the existence of corner singularities for solutions of the mod- 
ified Helmholtz equation in a polygonal domain. The basic differential form 
for the Modified Helmholtz equation is given by 



(22) W{z, z; A) = e-'^'+'l^''/^ 



dq . , \ J /i/? dq 



^^+iA, d.+ (-,--) dz- 



The associated global relation is given by W = 0. We concentrate on 
the particular corner Zi. Introducing the local coordinates {p,0) defined by 
z — Zi = pe^^, and using 

^ ^ dz~ 2 \dp pde)' dz~ 2 \dp^ pde 
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Figure 3. Gauge choice for comer singularity analysis. 



we find that 

(24) W = \f,~'^^''P^s)+ip^z{p,e)/\ 



Idq 

'~p~de 



+ ( \e'' + 



0^ 



dp. 



Without loss of generality we position the polygonal domain in the complex 
z-plane so that Zi = i, and \m Zj < for j ^ i as in Figure [5| Then, for z 
on any side other than on the sides Fj, Fj_i, and for A real and positive, we 
have 



(25) 



-iXz+il3^z/\ 



for some e > 0. Letting Zj = i, we find 

(26) e-i^^+'^'^"/^ = eA+/3VAg-pfc(A,e) j^^^^^-^ drf ^ g . 

For both Fj and Fj_i, ReA:(A,0) > 0, since 

-7r<6'<0, i.e. \e + 7r/2\ < tt/2. 
The sides Fj_i and Fj yield the following contribution to the global relation: 



(27) ie 



A+/32/A 



ir.l 



-pk{x-ei) 



Idq 



+ 



dp 



le 



|r.-i| 



Ae-'^'-i + 



Idq 

~p~de 

0" 



{P:-Oi-i)- 



j{p,- 



yi-ij 



dp 



We can examine the possibility of corner singularities by considering the 
limit of equation (27) as A — )• 00. I what follows we analyse several types of 
boundary conditions. 
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5.1. Continuous Neumann-Neumann. Suppose that Neumann bound- 
ary conditions p^'^dq/ dO{p, -de) for i = i — l,i are prescribed at the sides 
Ti-i and Fj. Furthermore, assume that the given data are continuous at the 
corner. Let 



Dep"^' as p ^ 0, 



where dg > —1 for ^ = i, i — 1, so that the integrals in (27) are weh-defined. 
Employing Watson's lemmsj^ in ( |27[ ) we find 

Ae-'^' AF(^^ + 1) ^ Ae-'^'-^ A-iF(d^-i + 1) 



(28) 

Hence 



i-l 



di, 



-lAi 



D 



Note that convexity implies that Aj^j_i = Oi — aj_i G (0,7r). We deduce 
that there exists an integer Ms{0,l,2,...}, such that 

2ttM 



(29) 



di-i — di 



D 



Di. 



i,i — 1 imply that M cannot be negative. 



1 are 0(1) as p | or 



The inequalities di > —1 for 
We deduce that either both g(/j, —6() for £ = i,i 

= di > 2. 

In either case, corner singularities do not occur, independent of the internal 
angle Ai,j_i. 

5.2. Continuous Dirichlet-Dirichlet. Suppose that Dirichlet boundary 
conditions q{p,—6i), £ = i,i — 1 are prescribed at the sides Fj and Fj_i. 
Furthermore, assume that the given data are continuous at the corner. Let 
dq 

de 

where rig > for 



-{p, -9e) ~ Nip 



as p 



0, 



I, I 



1, so the integrals in (27) are well-defined. 
The leading order term in A as A — )• oo originating from the Dirichlet data 
vanishes due to continuity. Indeed, this term gives rise to the contribution, 
as A ^ oo, 

i(q{Q,-ei)-q{Q,-ei^l)), 
which vanishes. Employing Watson's lemma in (27) we find 
NiT{ni) ^ iV,-iF(ni_i) 
(AeK-«'+'^/2))''^ (Aei(-^'-i+'^/2)) 

Thus, there exists an integer M G {1,2,...}, such that 

2-kM 



■rii-i • 



Ui-l 



rii 



A 



Ni^ 



Ni. 



i'i—l 



Employing similar reasoning as in the previous section we deduce that corner 
singularities do not arise, independent of the internal angle Aj^j_i. 



^Recall Watson's lemma e-'^p* dp ~ V{e. + 1) 
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5.3. Vanishing at the corner: Dirichlet-Neumann. Suppose that Neu- 
mann and Dirichlet boundary conditions are prescribed at the side Fj and 
Ti-i respectively. Furthermore, assume that the given data vanish at the 
corner. Let 

q{p, -Oi) ~ A/% |^(/0, -Oi-i) ~ Ni^ip""'-^ as p ^ oo, 



where di> —1 and nj_i > 0. Employing Watson's lemma in (27) we find 

\e-'^^DiT{di + 1) ^ iV,-iF(n,_i) 

Hence 

We deduce that there exists an M G {0,1,2,...} such that 

TT / _ 1 



Ni-i = Didi, di = m^i = 2M + ^ 

Singularities can occur only if nj„i < 1, i.e. if Aj^j_i > 7r/2. We conclude 
that corner singularities can occur only if there is an internal angle smaller 
larger than 7r/2, in agreement with the classical results for elliptic boundary 
value problems in corner domains [5]. 

5.4. Discontinuous Dirichlet-Dirichlet. In this case, as A — t- cxd, there 
exists a non-vanishing term proportional to the discontinuity. This implies 
the existence of non-integrable singularities. 

Remark 4. The limit as A — )• yields precisely the same information as 
the limit A — )• oo. This is not surprising, since (27) is invariant under the 
transformation A i— )• /3^/A followed by complex conjugation. 

6. Non-Integrable Singularities 

It was shown in ^that discontinuous Dirichlet boundary conditions can 
give rise to non-integrable singularities. In order to eliminate such a singu- 
larity at the corner Zj, we introduce the differential form V defined by 

V{p, A) = W{p, 0, A) + d (^ieA+/3VAe-pfc(A,^^) £ 0) dp' j , 

where 7^ 0. Along a ray of constant 0,V\s given by 

?2 \ rp 



(30) V{p, e, A) = ie^+^'/^e~^^(^ 



./ /32 \ fP Idq , , 



d/5 



and hence the non-integrable singularity is eliminated. A boundary value 
problem involving several non-integrable singularities, can be decomposed 
into a series of problems each involving a single corner. In what follows we 



20 



A.C.L. ASHTON AND A.S. FOKAS 



analyse in detail a particular boundary value problem which has two non- 
integrable singularities. In this case it is possible to eliminate both of these 
singularities using a single differential form. 

Example. Let q(x, y) satisfy the modified Helmholtz equation in the inte- 
rior of the half-strip formed by the corners 

zi = CO + i£, Z2 = ii, Z3 = 0, Z4 = 00, 

where i is real and positive. Let 

(31) q{x, 0) = q{x, £) = 0, < x < 00; q{0, y) = I, < y < £. 
The symmetry relation q{x, y) = q{x, i — y) implies 



dq 
dy 



(a:,0) = -|^(x,^), |^(x,^^)=0, < x < 00. 



It is possible to eliminate the non-integrable singularities at the corners Z2 

and Z3 using the following differential form: 

(32) 



V = W + d 
where 



0(A) 



dy 



x',y)dx' + e^^'-'+'^yc2 



dq 

dx 



X, y') dy' 



-ilX-^y u;(A) = A + ^ 



and ci, C2 are constants that satisfy ci — C2 = 1. Indeed, the definition of W 

dq 



in (22) implies 



dy 



+ ujq ] dx + 



dq 
dx 



nq dy 



Using this equation and employing the identities 



y dh. 



e/2 



dx"^ 



/ dq^ 



(x, y') dy' = 4/3^ / q{x, y') dy' - -^{x, y), 



1/2 



X q2 



d 1r / 



Qyi^.^' ^y)<^x' = ^P"^ I Q{x,y')dx' - ^{x,y), 



dy 
dx 



equation (32) becomes 



(33) W = e^^^'^y 



+ 



(rx Qq ry 
''Ldy^"'^'^'^"'^"'l/2d-x^"'''^'^'\ 
ry 

+ A(3'^C2 q{x,y')dy' +ujq{x,y) 
Je/2 

( /"^ dq fy dq 

.(^c,y^-(x',,)dx'+c2y^^^^(.,y')dy' 

/•X 

-F4/3^ci / q{x',y)dx' -Q.q{x,y) 
J 00 



dx 



dy 
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Thus, the Neumann boundary values, which involve the non-integrable sin- 
gularities, have been eliminated. 

For simplicity we set ci = 1 and C2 = 0, thus, (33) becomes 

(34) 



W 



— dx' + ujq{x,y) 



dx 



+ 



ijj 



dx' + 4/3^ / q{x' ,y)dx' -Q.q{x,y) 

' Joo 



dy 



dy 

Using this differential form and employing the method of [U [7] , it can be 
shown that q is given by 



(35) q{x,y) 



where 



1 



-flx-Loy G{X) dX 



+ 



-Qx—ujy 



G{X) 



dX 



^-nx+uj{£-y) G{X) dX 



l + e^'^ X 



G(A) = VL{X) 



In what follows we verify that the expression, as defined by (35), satisfies 
the given boundary conditions. 



Boundary condition q{x,Q) = 0.' Evaluating equation (35) at y = we find 
g(a;,0) 



1 



G(A)^ dA ^ 



l + e^'^ X 



e-^-G{X) ^ 

lOO ^ 



Noting that 



G{X) dX 



G{X) 



dX 



-nx+oj£ G{X) dX 
l + e"^^ X 

- G{X) dX 
1 + X 



we find 



q{x,0) 



1 

'2^ 





+ 

ioo J 



'Qx 



G{X) 



dX 
T 



0. 



Boundary condition q{x,i) = 0: This follows in a similar fashion. 



Boundary condition q{0,y) = 1: Evaluating (35) at a; = we find 



(36) g(0,2/) 



1 

'2tt 



-.y G{X) dX 
l + e^^' X 



-wy 



G(A) 



dA 

T 



+ 



G(A) dA 



l + e"^^ X 

In order to take care of the singularity at A = i/3, we deform the contour of 
the second integral on the right hand side of ( 36 ) to the contour L depicted 
in Figure 4. Then, q{0, y) can written in the form 
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i/3 



Figure 4. Contour L in the complex A-plane. 



(37) 



?(o,y) 



1 

'2^ 



+ 

L Jo 



uj A 



+ 

L Jo 

oo 



-ujy 



-2 I e 



-ujy_ 



OJ A 
J] dA 



1 + e^^ w A 

Jordan's lemma implies that the first integral on the right hand side of 



(37) vanishes. The integrand of the third integral remains invariant under 



the transformation A i— )■ 1/A, thus this integral also vanishes. The second 



integral on the right hand side of (37) has a pole at A = i/3 with residue — i. 
Hence g(0, y) = 1. 



The non-integrable corner singularity. Equation (35) implies 



1 



l\ dA 



e^^^ + lj A 

— oo 



+ 



„ A 



+ 



l\ dA 



10 ve'^^+iy ^ 

Using a similar argument to the one used in calculating q(x,0) we find that 
the expression on the right hand side simplifies to 



f 



f^(x',0) dx = - I e 
dy vr 7o 



-i(A-/3VA) ' ^ ' ' - i \ CLA 



3(A+/32/A)^ + 1 / A ■ 



We will show that 



(38) 



da 4 

7e(2;',0)dx' = logx + C'(l) asx|0. 

X dy TT 



Indeed, letting k = X — /3^/A, i.e. 

dA 1 



we find 
(39) 



^(x',0)dx' = - 
ay vr 



A + ^ = Vfc2 + 4/32, 
A 



-ifcz 



1 



+ 



ev- +1 

^ /•oo 



dA; 



V^A;2 + 4/32 
dA:. 
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The integrand in the first term is 0{l/k'^) as k — )• oo, so is absolutely 
integrable. Hence, the first term, being the Fourier transform of an element 
of L^(R), gives rise to a bounded, continuous function of x. The second 
integral is 

2 coskx „ 2 cosfe „ 2 f°° cosk 

- / dfc = - / = dfc + - / = dk. 

Jo ^/WTI^ ^ Jo x/A;2 + 4/32x2 vr ^A;2 + 4/32^2 

It is straightforward to show that the second of these integrals is a bounded, 
continuous function of x. Writing cos k = (cos — 1) + 1 we deduce 

-^ix',0)dx' = - +0 1. 

dy' ttJo VA;2 + 4/32x2 



The claim in (38) now follows by computing the elementary integral on the 
right hand side of the above equation. 



7. Appendix 

It was claimed in §4 that a distribution in iS'[0, 1] can be written as a 
finite sum of (weak)-derivatives of bounded functions supported on [0,1]. 
We have been unable to find the proof in the literature, so here we present 
a proof. 

Theorem. For each u £ f [0,1] there exist a finite collection of complex 
Borel measures {/xq, • • • , /^n} such that 

n „i 

{u,ip)=Y, ^^"'\T)dfim{T) V(^G^[0,1]. 

Proof. Let u G £'[0, 1] with ord(u) = n. Define the map 

r : £[0, 1] ^ ^[0, i]e- - -ec[o,i] (^^,^',..., , 

n + 1 copies 

where C[0, 1] denotes the Banach space of continuous functions on [0, 1] 
equipped with the supremum norm. Note that T(ip) is uniquely determined 
by if. Now consider the linear functional A : r(iS[0, 1]) — )• C defined by 

A{Ti^)) = {u,ip). 

This map is bounded when T{S[0, 1]) is treated as a subspace of C[0, 1] © 
• • • ® C[0, 1]. Indeed, since u £ £'[0, 1] has order n, we have 

iA(r((^))i = \{u,yp)\ <cj2 ii/'^iioo G m !]• 

By the Hahn-Banach theorem, A has an extension to all of C[0, 1] © • • • © 
C[0, 1], and by the Reisz-Markov theorem this extension must have the form 

n „i 

A{tpo,...,tpn)=y2 ^m(T")d/Xm(T"), 
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where {/iq, . . . , /in} are complex Borel measures. In particular, 



n »i 



which is the representation stated in the theorem. 



A similar result holds if one considers the map 



r:,f[0, 1] ^ L^[0,l]®■■■®L^[0,l]■ 



lI^ this case one finds a collection {/o, . . . , /„} in L°°[0, 1], such that 
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